Classical picture of post-exponential decay 
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Post-exponential decay of the probability density of a quantum particle leaving a trap can be 
reproduced accurately, except for interference oscillations at the transition to the post-exponential 
regime, by means of an ensemble of classical particles emitted with constant probability per unit 
time and the same half-life as the quantum system. The energy distribution of the ensemble is 
chosen to be identical to the quantum distribution, and the classical point source is located at the 
scattering length of the corresponding quantum system. A ID example is provided to illustrate the 
general argument. 

PACS numbers: 03.75.-b, 03.65.-w, 03.65. Nk 



I. INTRODUCTION 

Exponential decay is ubiquitous in quantum physics 
and constitutes the typical dynamical pattern for unsta- 
ble systems. Theory predicts deviations at short times, 
related to the Zeno effect, and also at long times, with sig- 
nificant implications for cosmology [l|, hidden variables 
non-hermitian formulations jj, 0]' radioactive- 
dating methods [5]. On the experimental side, the post- 
exponential regime, generally algebraic, has been elusive. 
There are few claims to have observed it @, which has 
engendered much effort to explain why, or to improve its 
observability 

Classically, exponential decay arises from a constant 
decay probability per unit time. What has been lacking 
so far has been a correspondingly simple, physically ap- 
pealing classical picture of post-exponential decay. The 
purpose of this short article is to provide such a frame- 
work, Iwfollowing up on an old suggestion of R. G. New- 
ton 

The standard quantum mechanical derivations of post- 
exponential decay are not very helpful in suggesting such 
a picture. The early reliance upon the Paley- Wiener the- 
orem jlTJ, is not very illuminating from a physical per- 
spective and it provides bounds, not precise predictions. 
Similarly post-exponential decay can be attributed math- 
ematically to the fact that the pole contribution to ex- 
ponential decay is eventually comparable to or smaller 
than a line integral, whose value arises predominantly 
from a saddle point at threshold, associated with slow 
particles 0,[i3,[il|- This is surely more intuitive, yet not 
fully satisfying for those seeking a pictorial, rather than a 
complex variable, understanding of the phenomenon jl^ . 
In this vein, Hellund proposed an electrostatic analog 
[Is'] which relates quantum emission of radiation to the 
damped oscillation of a charge describable in purely clas- 
sical (stochastic) terms, and interprets the deviation from 
exponential decay as a "straggling phenomenon" charac- 
teristic of a diffusion process. However, quantum dynam- 
ics cannot generally be reduced to a classical diffusion 
process. Jacob and Sach fld\, in their field-theoretical 



analysis of a scalar particle coupled to two pions, found 
nonexponential terms decaying like t~^^^ in the ampli- 
tude. Their explanation was geometrical: a particle pro- 
duced at position r having velocity between v and dv, will 
appear after a time t within a spherical shell of radius vt 
centered on r, and thickness tdv. The probability that it 
will be found in a small volume element within the shell is 
inversely proportional to its volume Ant'^v'^dv. Hence the 
probability amplitude is proportional to t~^/'^. This can- 
not be a universal explanation though, since the decay 
amplitude in ID models behaves generically like t~^/'^, 
whereas the above argument translated to ID would im- 
ply only 

Post-exponential power-law behavior is sometimes in- 
terpreted as expressing the dominance of free-motion 
[itI [isj . However explicit calculations of the long-time 
propagator for specific potential scattering models show 
that in general both free and scattering terms are needed 
in the propagator, to reproduce the correct results jl9l] . 

A general argument that justifies non-exponential de- 
cay is the so-called "initial state reconstruction" [^O, [2l| . 
Deviations from exponential decay would be the conse- 
quence of Feynman paths returning to the initial state 
from the orthogonal, decay products subspace. This ar- 
gument alone, however, makes no quantitative predic- 
tions of the observed behavior and is applied to the sur- 
vival probability, which is not always easy to observe. 

In this paper we take up and develop a frequently over- 
looked observation by R. G. Newton making use of classi- 
cal mechanics Newton noted that if a point source 
emits classical particles with an exponential decay law 
and with a suitable velocity distribution, the current den- 
sity away from the source will eventually depend on time 
according to an inverse power law. Indeed, we shall show 
that by adjusting the parameters according to the quan- 
tum system, the classical model accurately reproduces 
the onset, power law, and intensity of post-exponential 
decay of the quantum probability density of a particle 
escaping from a trap. For simplicity we assume that the 
particle is restricted to the half-line, r > 0, as in quantum 
s-wave scattering. We shall also assume that the initial 
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quantum state is orthogonal to any bound states so that 
it must eventually decay (escape) fully from the trap. 



II. CLASSICAL AND QUANTUM SOURCES 
AND DECAY 



Consider first a source at ?■ = which emits classical 
particles with a definite velocity v from to = 0, so that 
the fraction of particles emitted between to and to + dto 
is 



Pito)dto - ^e-*°/^ 



(1) 



where r is the emission life time. The spatial probability 
density observed at point r, at time t, is 



Pc 
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(2) 



where 9 is the step function. For the more general case in 
which the emitted particles have a velocity distribution 

Pcir,t)= dvp{v)—e-^'--l-^'\ (3) 

Jr/t VT 



or, using r = vit — to), 



-ta/T 



(4) 



The behavior for t >>> t of this integral can be ex- 
pressed as an asymptotic series, 

m 

P,{r,t) ^ ^r"[g(")(0) -g(")(t)e-*/l, (5) 



n=0 



where 



1 / r 
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(6) 



and gr^") is its n-th derivative with respect to to- The 
leading term asymptotically is 



Pr 



(7) 



This is equivalent to Newton's result [§, [13] (we use the 
probability density rather than the current density). To 
advance from here, consider now the decay from a quan- 
tum trap of a system prepared in a normalized non- 
stationary state I'^o)- The wave function of this state 
at a point r and time t is 



^ir,t)^{r\e-^'''/^\^o), 



(8) 



with corresponding probability density Pq (r, t) ~ 
4'(r, t)|^. Using stationary states normalized in energy 



usij) (such that {u^'\ue) = 5{E — E )), and inserting 
the completeness relation, 

/•OO 

^{r,t)^ dE{r\uE){uE\^o)e-'^''''. (9) 



The ue are solutions of the s-wave, radial Schrodinger 
equation. 



UE{r) = 0, 



(10) 



where v{r) = {2ni/h^)V{r) and fc^ = {2m/h^)E. As 
in 22], it is convenient to define new solutions Wk{r) ~ 

h'fJ^UEir) normalized as (wj./ ]i/7fe) — S{k — k), which 
obey the boundary condition 



lim Wk{r) = y — sin[fcr- 



S{k)], 



(11) 



where 6{k) is the phase shift of the s-partial wave. These 
solutions are related to the regular solutions 0^ (r) (which 
behave like the Riccati-Bessel function jo{kr) as r 0), 
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(12) 



where f{k) = \f{k)\exp{~i6) is the Jost function, as 
defined for example in Taylor |23j]. It gives the rel- 
ative normalization between solutions having unit in- 
coming flux at infinity, and solutions that have slope 
k at the origin. The partial-wave S'-matrix element is 
S{k) — f{—k)/f{k). Zeroes of f{k) in the upper half 
complex momentum plane correspond to bound states, 
while those in the lower half plane are associated with 
scattering resonances. 

For an initially localized non-stationary state 



(rj^'o) =0 for r > Ta, 
the wave function can be written as 



*(r,t) 



2m 



iEt/h 



(13) 



(14) 



which has a form similar to the survival amplitude ob- 
tained in . They have generically the same asymptotic 
behavior at long times, which corresponds to an energy 
distribution p{E) = \{ue\'^o)\'^ E^/^, as E ^ 0. This 
long time asymptotic behaviour is governed by the prop- 
erties when k oi the integrand of Eq. (fT4|). For 
"well behaved" potentials, (those falling off faster than 
when r — )> oo and less singular than r~^^^ at the 
origin) , the £ — Jost function tends to a constant when 
fc ^> 0. (In the exceptional case that a zero energy res- 
onance occurs f{k = 0) = 0.) The 0^ behave near the 
origin as Ricatti-Bessel functions, jo{kr), and are there- 
fore linear in k. The behaviour of the integrand near 
threshold is thus ^ E^^^. Following the same steps as in 
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the derivation of the asymptotic behavior of the survival 
amplitude in |22|] , one finds that the position probability 
density behaves like Pq{r,t) ^ t^^ at long times. 

The energy distribution corresponds asymptotically to 
a velocity distribution since all particles are eventually 
released. The two distributions are related by 



g{E)dE = piv)dv. 



(15) 



Setting E = ^mv'^, m being the mass of the emitted par- 



ticles, makes g{E) ^ E^l"^ =^> p[v) 



Going back to 



Eq. ([7]) and considering the long time regime t >>> r, 
the classical particle velocity can be approximated by 
V = r/t, so p{r/t) = p{v), which implies, as expected, 
that at large t the main contribution to the position prob- 
ability density is from slow particles. If we consider the 
same dependence as in the quantum case, p{v) ^ w^, Eq. 
([7]) implies an asymptotic behavior Pc{r,t) ~ t^^, i.e., 
the classical model leads to the same power law depen- 
dence as the quantum one. Moreover, in the following we 
shall see that it can be adjusted to provide the correct 
amplitude factor as well. 

Let us return to Eq. (ITU) and write the bra-ket factor 

as 



(16) 



ik\^o) = \i^\f{mwk\^o), 



Using Eq. (|12p and the asymptotic behavior given by Eq. 
(|lip. the wave function for r — > oo may be written as 



*(r,t) 



sinffcr 



(5(fc)]e~*^*/''. (17) 

At low energy, the phase shift 6{k) is well described by 
the effective range expansion 

fccot(5(fc) = — - + -rok'^^ , (18) 

where oq is the scattering length while vq is called the 
effective range of the potential function. The asymp- 
totic form of the resulting integral can be obtained from 



the Riemann-Lebesgue lemma 2j|. Only the main term 
which depends on the fc — s- behavior of the integrand is 
kept. This gives a probability density of the form 



1 



(19) 



where /3 is the strength factor for the asymptotic depen- 
dence of the velocity distribution, p{v) which will 
depend on the particular state and potential. It has units 
of [/3] ^ . In the approximation v — r/t we can write 



p(.)^/?-. 



(20) 



Introducing this in Eq. ([7]) we obtain for the classical 
probability density at long times 




\&)-y„g(r-L) 




FIG. 1; (Color online) Scheme of the Winter model, a) The 
initial state is the ground state of an infinite square well, b) 
One of the walls is replaced by a delta barrier V — Vo5{r~L). 



Compare now Eqs. (IT^ and (PT|) . and to avoid confusion, 
let us rewrite r ^ rq for the quantum case and r ^ for 
the classical one. We see that if the classical coordinate 
is shifted by oq, r,. — Vq — ao, the classical model will re- 
produce the quantum probability density. Equivalently, 
Pc{r) — Pq{r) at long times if the classical source is not 
at the origin but displaced by the scattering length qq. 

In the exceptional case of a potential with a zero energy 
resonance ao — > oo, and therefore the first term on the 
r.h.s. of Eq. ([18]) is absent. This causes the Jost function 
to have a simple zero at fc = 0, see 23], and therefore 
|(wa;=o|^'o)P is nonvanishing. We then find, instead of 
Eq. that 



Pq{r,t) ^ \{wk=o\'fo)\^m/{ht), 



(22) 



which is again in agreement with the classical expression 
dZl) taking piv = 0) = \ {wk=o\'^o)\^m/h. 



III. MODEL CALCULATION 

Now we check the above general results for Winter's 
delta-barrier model [l^ which is described in Fig. [T] 
The initial state is an eigenstate of the infinite square 
well potential. 



(^l*o) 



Esin(^) r<L 



r>L 



(23) 



and the /c-normalized basis functions are 



{r\wk) = e-''^(^) 



i{kr)/f{k) 



r < L 



(V2)[e- 



-ikr 



S{k)e'''''] r>L' 
(24) 

where S{k) = f{k)*/f{k), and the Jost function for this 
model is 



(21) 



(25) 
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with a = 2mVb/fi.^. p{v) = \{wk\^o)\'^m/H can be calcu- 
lated exactly and takes the form 

Lm 

p{v) 



nh fc2 + ak sin(2fcL) + sin^(fcL) 
2m: sin(fcL) 



fc2L2 



(26) 



From here the exact classical probability density is cal- 
culated numerically using Eq. whereas the quantum 
density is given by the square modulus of Eq. In 
the large-i region the probability density has analytical 
expressions in both quantum and classical cases given by 
Eqs. (IT^ and (PT|) respectively. The coefficient (3 is easy 
to find from Eq. (l26l) in the limit u — t- 0, 



f3 



3 r3 



(27) 



(1 + aL)^nVh^' 

Also, Eq. dm) and S{k) = e^^^^''\ give for the Winter 
model the explicit source shift 

aL2 



ao 



1 + aL' 



(28) 



which, for a > 0, lies between (for aL 0, no barrier), 
and L (for large aL, strong confinement). 

Finally, the quantum and classical probability densities 
both take (shifting the classical point source by ap) the 
post-exponential form 



Lm\ ( aL^ \ 1 

n2(l-)-aL)2 \~^J \^ l + aLj 

(29) 

The agreement is illustrated in Fig. [2J where the exact 
decay curves (numerically integrated) are plotted. The 
classical density (triangles) indeed reproduces the quan- 
tum behavior (solid line) if the source shift is taken into 
account. For comparison we also show a curve in which 
the shift is not applied, so that the classical source re- 
mains at r = (circles). Taking the same value for r, we 
see that the classical model also agrees with the quan- 
tum one in the pre-exponential and exponential zones 
(0 < t < 5 in the drawing). The classical model dif- 
fers only in the absence of oscillations which occur at the 
onset of post-exponential behavior, due to quantum in- 
terference. The asymptotic behavior is indistinguishable 
on the scale of the figure from the analytical expression 
Eq. (EID. 

The exceptional case of a zero energy resonance cor- 
responds to an attractive delta with a — —1/ L. In this 
case from Eqs. (|7l22p we get 

4Lm 



(30) 



for both the classical and quantum cases: see Fig. [3l If a 
is close to the critical value, say a — —1/L + e, the decay 
follows a decay law for some substantial period of 
time until the decay eventually dominates, see Fig. 
m The smaller is e, the longer does the behaviour 
persist. 




10 20 30 40 50 
t 



FIG. 2: (Color online) The exact probability density obtained 
numerically versus t: quantum result solid (red) line; classi- 
cal solution including the source shift, Eq. ([28} (triangles) ; 
classical solution without applying the source shift (circles). 
Parameter values: h = 1, m = 1/2, L — 1, a = 5, n = 1, 
r = 2, r = 0.5. 
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FIG. 3: (Color online) Logarithm of the probability density 
versus t, for an attractive delta potential having a zero energy 
resonance: exact quantum numerical solution (solid red line); 
classical model solution (triangles). The long-time behaviour 
is indistinguishable from Eq. (|30|l . a — —1, r — 0.2, and 
other parameters as in Fig. 2. 



IV. DISCUSSION 

To summarize, the above results provide an intuitive 
physical picture and quantitative description of post- 
exponential decay of the probability density at points 
distant from the source. We have developed the classical 
model suggested by Newton so as to achieve an accurate 
match between classical and quantum decays. Purely ex- 
ponential decay from a source leads naturally, because of 
dispersion associated with a velocity distribution of the 
emitted particles, to the same power law decay in quan- 
tum and classical scenarios. Quantum mechanics is re- 
quired to provide the emission characteristics, but initial 
state reconstruction (ISR) plays no role in the classical. 
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FIG. 4: (Color online) Logarithm of the probability density 
versus t for a delta potential slightly less attractive than re- 
quired to produce a zero energy resonance: exact quantum 
numerical solution (solid blue line), and approximate power 
law decays proportional to (Eq. (|30p . dashed red line) 
and f"^ (Eq. dotted green line), a = -0.98, r = 10, 

r = 0.2, and other parameters as in Fig. 2. 



r < L ) is still in its exponential regime when the tran- 
sition shown in Fig. 2 (at r = 2) takes place, i.e., the 
purely exponential decay hypothesis ([T]) for the classical 
source is justified, and the onset of the post-exponential 
regime of survival within the trap cannot causally affect 
the transition observed in the density outside the source. 

Due to recent advances in lasers, semiconductors, 
nanoscience, and cold atoms, microscopic interactions 
are now relatively easy to manipulate, decay parameters 
have become controllable, and post-exponential decay 
more accessible to experimental scrutiny and/or applica- 
tions 0. Under appropriate conditions it could become 
the dominant regime and be used to speed-up decay via 
an Anti-Zeno effect ^26^]. Moreover, recent experiments 
on periodic waveguide arrays provide a classical, electric 
field an alog of a quantum system with exponential de- 
cay [13, l28l | where the post-exponential region could be 
studied in a particularly direct way. 
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